Inflation in Oscillating Universe 
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Abstract 

We make use of possible high energy correction to the Friedmann equation to implement the 
bounce and study the behavior of massive scalar field before and after bounce semianalytically 
and numerically. We find that the slow-roll inflation can be preceded by the kinetic dominated 
contraction. During this process, the field initially in the bottom of its potential can be driven by 
the anti-frictional force resulted from the contraction and roll up its potential hill, and when it 
rolls down after the bounce, it can driven a period inflation. The required e-folds number during 
the inflation limits the energy scale of bounce. Further unlike that expected, the field during the 
contraction can not be driven to arbitrary large value, even though the bounce occurs at Planck 
scale. 



PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 



The idea of oscillating universe is ancient, in which the universe oscillates through many 
successive expansion/contraction cycles, which was proposed originally by Tolman l| in the 

n 

1930's. Recently this idea has been awaken again in cyclic scenario jyj, which is motivated by 
the string/M theory. The relevant dynamics with primordial perturbations can be described 
by an effective theory in which the separation of the branes in the extra dimensions is 
modeled as a scalar field. But the problems around bounce may remain. For contracting 
universe, it is hard to escape from a singularity in the frame of general relativity. Thus for 
the scenario of oscillating universe, one of difficulties of implementing it is how to obtain 
a non singular bounce. Several proposals have been discussed, such as using the negative 
energy density ODD or curvature term p, or some high order correction term in the 
action [7, 0]. 

During contraction of each cycles, the scalar field can be driven and roll up along its 
potential to large enough value for inflation to occur, which has been noticed in Ref. js( for 
closed universe, in which it has been shown that an oscillating universe certainly undergoes 
inflation after a finite number of cycles. During this process, a contracting phase dominated 
by kinetic energy followed by an inflation. Further, in Ref. ^(|, such scenario has been 
proposed to provide a possible explanation for the observed low CMB anisotropics on large 
angular scale. In this proposal, a kinetic dominated contracting phase is matched to an 
inflation phase. The power spectrum during the contraction ~ k 3 , which lead to an intense 
suppression of CMB quadrupole, in the meantime the nearly scale-invariant spectrum from 
slow-roll inflation is recover on small scale. But a physical mechanism for bounce is not 
included. In this paper, we make use of possible high energy correction to the Friedmann 
Eq. |l5j] to implement the bounce. We firstly focus the behavior of massive scalar field 
during such a bounce semianalytically and find that the slow-roll inflation can be preceded 
by the kinetic dominated contraction and the required e-folds number during the inflation 
after the bounce limits the energy scale of bounce. Finally we check and confirm these 
eatures numerically. Note that a similar scenario has been studied in string inspired frame 
1^ | . in which the Gauss-Bonnet term is used to construct a non singular bounce and the 



1 This is also the usual results of Pre Big Bang scenario 0] , see Ref. ^3 f° r 



a review. 
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superinflation regime before slow-roll inflation leads to a suppression of primordial spectrum. 
This has also a similarity to the case in loop quantum gravity jl^ . The effects from loop 
quantum gravity can driven the inflaton to its potential hill and increase the parameter space 



latory universe induced 
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of initial conditions for successful inflation. Furthermore, the osci 
by the effects from loop quantum gravity was investigated in Ref. 

II. DYNAMICS OF MODEL 

In this section we study the cosmological behavior of massive scalar field before and after 
a bounce phase semianalytically. We use the modified Friedmann equation 

h 2 = \{p v - P h (1) 

to implement a realistic bounce, where 8ir /m 2 = 1 has been set, p v is the energy density 
of scalar field and a is the bounce scale, when p v = a, the universe bounces. Eq. (JTJ) may 
be motivated in brane world scenario 1 511 . see also [17|. where our universe is embedded 
in high dimension space/time, see Ref. Il8|l for a recent introduction. Further it has been 
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pointed out in Ref. that due to the effects of the bulk and Israel matching conditions 
there might be some forms h(p) which is not standard Fridmann-like. 
The motion equation of scalar field <p with the mass m is 

(p + 3h(p + m 2 <p = (2) 

We assume that initially the universe is in contracting phase and the field ip is in the 
bottom of its potential. Thus p^ is small and the term p 2 , in Eq. is negligible. Some 
possible fluctuations will make (p departure from its minimum and oscillate near tp = 0. Since 
in this case m 2 3> h 2 , i.e. the frequency of oscillation is much larger than the evolution rate 
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of universe, thus after implementing the change of variable ip = a~2u, 



u + m 2 u ~ (3) 
can approximately be obtained, which can be solved as an oscillation 

u ~ u a sin {mt) (4) 
where u a is a constant amplitude. Thus we have 
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cp = a~zu a sin (mt) (5) 



When taking the time average over many oscillations of field, < >~ is obtained, which 
is similar to the case dominated by matter. Thus in this regime the universe contracts as 
a ~ if: 2 / 3 and the energy density of the field grows as 



p v ~ a 3 u a m (6) 

From (jUJ), we can see that with the contraction of universe the energy density of scalar field 
increases. Thus when h 2 ~ p v > m 2 , the term relevant with h of Eq. (J2J) can no more be 
neglected. In this case, a~ 3 u 2 m 2 > m 2 , which can be reduced to ip 2 a > a~ z u 2 a > 1. Thus as 
long as the oscillating amplitude of ip field is larger than 1, the oscillation of field will end 
and the field will roll up along its potential. Instead of Eq. (jSJ), 

yi + 3!vp~ 0. (7) 

is satisfied. From Eq. ((7j), we have 

(8) 

where c is the integral constant. Thus with the further contraction the universe will enter 
into the regime dominated by kinetic energy of field rapidly, i.e. (ft 2 ^> m 2 (p 2 , where p v ~ p v . 
This leads to p v ~ Lp 2 ~ 1/a 6 , thus we have 

a 3 -cJ^(t s -t) (9) 



2 , / t s tfcs 



where t s is the the integral constant, and 

^^ + Vi 1,1( t^r ) (10) 

where (pk and tk are the value of field and time just entering into the regime dominated by 
kinetic energy, respectively. 

When the kinetic energy of if approaches a, we have to includes the p 2 correction in Eq. 
(HJ). Thus instituting Eq. (JHJ) into Eq. (JTJ), we obtain 

whose solution can be obtained as follows 



a 
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3c 2 



c 
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(t - h) 2 + — (12) 



2 y ' 2a 



where tb is the integral constant determined by a as well as the value of (p and a at t = tk- 
When t = tb, the scale factor of the universe arrives at its minimum and after tb the universe 
bounces and expands. 

For the universe dominated by stiff fluid with the same state equation p = p as the case 
dominated by kinetic energy of scalar field, its evolution will symmetric before and after the 
bounce, and up to all time. But the evolution of the universe driven by massive scalar field 
is different from that with fluid. After the bounce, the field <p will still roll up. However, 
since h > 0, the term 3h(p of Eq. (J2J) serves as a damping term. Thus the roll- up motion of 
ip will decay quickly. When the velocity of ip arrives at 0, it reverses and rolls down along 
the potential. For the case of ip > 1, we have m 2 < h 2 . The field will enters the slow-roll 
regime in which the universe is dominated by the potential energy of the scalar field and 



is satisfied. In this case p v ~ — p v and the universe inflates. The motion of field ip is given 



where ip r is the value of field at reverse. When m 2 ~ h 2 , i.e. <p ~ 1, the inflation ends and 
the field (p will reenter into an oscillatory stage. 

Therefore from above discussion, we can see that before the bounce, the field initially in 
the bottom of its potential can be driven by roll up its potential hill and after the bounce it 
reverses and rolls down, and drives the inflation of universe. 

III. NUMERICAL RESULTS 

Then we study the above evolutive process numerically. We take m = 0.001 and a = 
0.0001, and a(0) = 1, f{0) = and <p{0) = 0.0001 for initial values. The evolutions of a 
and h are plotted in Fig. 1. We can see distinctly that the bounce is implemented. In 
Fig. 2, we plot the evolution of field (p. Initially the field oscillates about the bottom of its 
potential many times, then the field ip roll up along its potential rapidly. The field still rolls 
up after the bounce, but its velocity will approaches quickly. Thus shortly it will reverse 
and roll down, and then enter into slow-roll regime. When the field decreases to ~ 1 linearly, 
the slow-roll ends and the field will oscillates about <p — again. These numerical results 



3h(p + m 2 (p ~ 0, 



(13) 



by 




(14) 
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are consistent with our semianalytical discussions in last section. The modification scale 
a determines the energy scale of bounce, thus in some sense also determines the maximal 
value (p r which the ip can roll up to, see Fig. 3 for numerical solutions with various a. We 
can simply estimate the value of ip r . When the field just enter into the stage dominated by 
kinetic energy, <jp\ ~ vn?Lp\ can be satisfied approximately, and when the field is at the point 
of the bounce, its value is (p\ ~ a. Thus from Eqs. (JBJ) and 0, we have 



/\2 /\6 / • \ 2 02 

t s -t b \ ( a b \ ( cp k \ m z Lp% 



(15) 



Vr<Vk + \l-ln(-—^) (16) 



t s — tie J \a k J \(p b J a 
Instituting it into Eq. (fTUjl. 

7T ln (— f a , 
3 

can be obtained, where the second term has been doubled since during kinetic dominated 
regime the evolution of the field can be regarded approximately as symmetrical before and 
after the bounce, which can also be seen from Fig. 4. For a reasonable estimation, we take 
(fik ~ 1, thus obtain <p T < 4 for a = 0.0001 and (p r < 10 for a = 0.1, which is compatible 
with the numerical results of Fig. 3. In this case there is a up-limit <p r < 11 for a — 1, 
in which the modification of Friedmann equation is at Planck scale. Therefore, from Eq. 
|T6j) . for massive scalar field initially in the bottom of its potential with a small velocity 
from fluctuations, the maximal value that it can roll up to during the contraction is only 
dependent on its mass m and the modification scale a. The larger (p r is, the longer the time 
that the universe after the bounce inflates is. 

Focusing on the chaotic inflation model j^], instead of = 1 on whole paper, we 

take rrip = 1 here. The e-folds number is given by 

N ~ 2tt^ (17) 

and Eq. (|T6*j) is modified as 



<Pr <^ + \/ T ^-ln(- 2 ^ T ) (18) 
V 127T m z (pi 

For iV > 60, this requires ip r > 3. Considering m ~ 10~ 6 required by observational primor- 
dial perturbations, the limit to a is given by o > 10~ 8 . Further for o = 1 the maximal value 
that the field can roll up to is given by <p r < 1 + 12 In 10 ~ 5.5, which is suitable for a 
successful inflationary cosmology but is far away from eternal inflation 22]. 
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FIG. 1: The scale factor In a (solid line) and the Hubble parameter h (dashing line) as a function 
of time. The inset is a zoom before and after the bounce, h has been multiplied by 1000. 




FIG. 2: The scale factor In a (solid line) and the field ^asa function of time. 




FIG. 3: The field 99 as a function of time for various a. From long dashing line to solid line, a 
takes 0.0001, 0.001, 0.01, 0.1, respectively. Since m is taken same, initially the oscillation of fields 
has same frequency, which can be seen from Eq. 
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FIG. 4: The potential energy (solid line) and the kinetic energy (dashing line) of the field as a 
function of time before and after the bounce. 




FIG. 5: The scale factor In a (solid line) and the field tp (dashing line) as a function of time for an 
oscillating universe. 
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FIG. 6: The potential energy (solid line) and the kinetic energy (dashing line) of the field as a 
function of time for an oscillating universe. The time taken only includes two cycles. 
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FIG. 7: The scale factor In a as a function of time for an oscillating universe. The long, short 
dashing and solid lines are the cases of = 50,5,2.5, respectively. For ^ = 2.5, from (|16l) . we 
have cp r < 1. In this case, inflation dose not occur in each cycle, thus the amplitude of successive 
cycle dose not increase. 

We have found that the slow-roll inflation after the bounce can be preceded by the kinetic 
dominated contraction. Further we would like to check this feature in a controlled model of 
oscillating universe 2 . We plot the evolutions of (a, (p) and the (kinetic, potential) energy 
with time in Fig. 5 and Fig. 6, respectively. To show the oscillation of universe better, 
we take m 2 = 0.2, a = 50 and a(0) = 1, p(0) = an d (p — 1 for initial values. We see 
that the maximum of expansion amplitude in each successive expansion/contraction cycle 
grows gradually. This asymmetry of cosmological evolution seems conflicted with intuition, 
since both the field equation (0) and Fridmann equation are non dissrpative and time 
reverse invariant. The reason of asymmetry has been analyzed in Ref. J9( for the case of 
closed universe. However, in fact this is not dependent on the curvature of universe, but 
the evolutive behavior of scalar field. For ip > 1, during the expansion, the state equation 
p v ~ — p v , thus the field is in slow-roll regime and the universe expands exponentially i.e. 
a ~ e ht . However, during the contraction, the universe is dominated by kinetic energy of 
scalar field. In this case, a ~ t 1 ^ 3 . Thus the time when the field rolls up to the value from 
which it will roll down into slow-roll inflationary regime is more shorter. It is the evolutive 
asymmetry of the field that leads to the growth of successive expansion/contraction cycles. 

2 To realise a transition from expansion to contraction, we take the density of potential energy < in the 
minimum. The details of cosmology with negative potential have been studied in Ref. [j^f . 
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When that m 2 is the same order as a, from ([160. we obtain </? r < 1. In this case, the field 
can not roll up to slow-roll regime. Thus the asymmetry between the expansion driven by 
the potential of field and the contraction driven by the kinetic term of field will not exist. 
We can see from Fig. 7 that the amplitude of successive cycle does not increase in this case. 

IV. CONCLUSION 

We make use of possible high energy correction to the Friedmann equation to implement 
a bounce. We study the behavior of massive field and find that the slow-roll inflation can 
be preceded by the kinetic dominated contraction. During this process, the field initially 
in the bottom of its potential can be driven by the anti-frictional force resulted from the 
contraction and roll up its potential hill. We show that the required e-folds number during 
the inflation after the bounce limits the energy scale of bounce, for iV > 60, o > 10~ 8 is 
required. Further unlike that expected, in our case the field can not arrive arbitrary large 
value, even though the bounce occurs at Planck scale, for a — 1, we have (p r ~ 5.5. The 
reason is that during the contraction the kinetic energy of the field is much larger than its 
potential energy and the growth of the field value is only logarithmical. When the kinetic 
energy approaches a, its potential energy is still far less, see Fig. 3. These features have 
been checked and confirmed numerically. Moreover, for possible oscillating universe model, 
we point out that the successive increasing of expansion amplitude of each cycle is actually 
depended on whether there exists inflation in each cycle. For the oscillating universe without 
inflation in its each cycle, its amplitude dose not increase with time. 
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